Draft version January 14, 2010 

Preprint typeset using I4TgX style emulateapj v. 1 1/10/09 



O 
G 



o 

u 

43 

Oh- 

6 

CO 



> 
in 

(N 

o 
o 



13 



THE MYSTERIOUS MERGER OF NGC 6868 AND NGC 6861 IN THE TELESCOPIUM GROUP 

M. E. Machacek, E. O'Sullivan, S.W. Randall, C. Jones, and W. R. Forman 

Harvard-Smithsonian Center for Astrophysics 
60 Garden Street, Cambridge, MA 02138 USA 
Draft version January 14, 2010 

ABSTRACT 

We use Chandra X-ray observations of the hot gas in and around NGC 6868 and NGC 6861 in the Tele- 
scopium galaxy group (AS0851) to probe the interaction history between these galaxies. Mean surface bright- 
ness profiles for NGC 6868 and NGC 6861 are each well described by double /3-models, suggesting that they 
are each the dominant galaxy in a galaxy subgroup about to merge. Surface brightness and temperature maps 
of the brightest group galaxy NGC 6868 show a cold front edge ~ 23kpc to the north, and a cool 0.62 keV 
spiral-shaped tail to the south. Analysis of the temperature and density across the cold front constrains the 
relative motion between NGC 6868 and the ambient group gas to be at most transonic; while the spiral mor- 
phology of the tail strongly suggests that the cold front edge and tail are the result of gas sloshing due to the 
subgroup merger. The cooler central region of NGC 6861 is surrounded by a sheath of hot gas to the east 
and hot, bifurcated tails of X-ray emission to the west and northwest. We discuss supersonic infall of the 
NGC 6861 subroup, sloshing from the NGC 6868 and NGC 6861 subgroup merger, and AGN heating as possi- 
ble explanations for these features, and discuss possible scenarios that may contribute to the order of magnitude 
discrepancy between the Margorrian and black hole mass - a predictions for its central black hole. 
Subject headings: galaxies: clusters: general - galaxies:individual (NGC 6868, NGC 6861) - galaxies: inter- 
galactic medium - X-rays: galaxies 



1. INTRODUCTION 

One of the most important questions facing models of 
galaxy evolution today is how central supermassive black 
holes, found at the centers of most galaxies, co-evolve with 
their host galaxies. The observed correlations between cen- 
tral black hole mass and stellar velocity dispersion in present 
epoch galaxies (Mbh-c*, Gebhardt et al. 2000; Ferrarese 
& Merrit 2000; Tremaine et al. 2002) and between black 
hole mass and galaxy bulge mass or bulge luminosity (Mbh - 
M b uig e ; Magorrian et al. 1998; Hitring & Rix 2004) provide 
striking evidence for this coevolution. However, the shapes 
of these correlations, particularly at the extremes of velocity 
dispersion or luminosity, are not well known, and, which, if 
either, of these correlations is a reliable predictor of black hole 
mass in galaxies with high central stellar velocity dispersion 
or large luminosity remains controversial (Lauer et al. 2007). 
Key to resolving this controversy is understanding the dynam- 
ical connections between galaxy interactions, the feedback 
cycle from active galactic nuclei (AGNs), and black hole fu- 
eling and growth. 

Recent observations show that galaxy evolution occurs pre- 
dominantly in moderately massive galaxy groups at redshifts 
> 1 (Cooper et al. 2006). Thus to build dynamical mod- 
els for the transformation of galaxies and the growth of their 
central black holes, we must understand how galaxy interac- 
tions and feedback affect gas in and surrounding galaxies in 
the group environment. Questions that remain poorly under- 
stood include: What dynamical processes control the mass 
and energy flows into or out of the merging galaxies either 
fueling the central AGN and/or initiating star formation, or 
quenching star formation and heating and enriching the in- 
tragroup medium (IGM)? When and how are the dark mat- 
ter and hot gas halos of the merging galaxy separated from 
the galaxy to become part of the group gravitational poten- 
tial? How efficiently are the outer stellar halos of interacting 
galaxies stripped? Are central stellar velocity dispersions of 



the galaxies affected by these interactions and on what dy- 
namical timescales? What role do non-hydrostatic group gas 
(IGM) motions, induced by galaxies and galaxy subgroups 
passing through the group core, play in the co-evolution of 
the system? 

Merging galaxies in nearby galaxy groups are particularly 
important laboratories in which to address these questions, 
since only in these groups do X-ray observations have suf- 
ficient angular resolution to directly observe the X-ray edges, 
outflows, cavities, shocks, ripples and tails characteristic of 
the dynamical processes at work on the galaxies and within 
the group. Measurements of temperatures and densities in 
these features allow us to constrain the velocities, orbits and 
interaction history of the galaxies, as well as probe the cycle 
of nuclear activity that transports matter and energy between 
the galaxy's nucleus and the surrounding gas. Modeling the 
effects of these interactions on nearby galaxies, that lie on the 
extremes of the black hole mass scaling relations, provides a 
unique opportunity to understand the galactic processes that 
determine the shape and validity of these relations. 

In this paper we present the results of Chandra observa- 
tions of the two brightest elliptical galaxies, NGC 6868 and 
NGC 6861, in the Telescopium group (AS0851). AS0851, 
shown in Figure [1] is an Abell class galaxy group 28 Mpc 
distant (Tonry et al. 2001), containing > 11 member galaxies 
(Garcia 1993). Both the brightest group member, elliptical 
galaxy NGC 6868 (20 /, 09 m 54.1 s ,-48°22'46"4, B T = 11.66, 
NED) and the second brightest group galaxy, SAO galaxy 
NGC 6861 (20 ,! 07 m 19.5V48 o 22'12:'8,B r = 12.12, NED), 
25 '7 to the west, show signs of recent interactions and are 
likely to merge. Surprisingly, Wegner et al. (2003) mea- 
sured a central stellar velocity dispersion of the less luminous 
galaxy (NGC 6861) of ~ 414km s" 1 , a factor 1 .6 greater than 
the 250 km s" 1 central stellar velocity dispersion of the more 
luminous elliptical galaxy NGC 6868. 

We summarize the galaxy properties and estimates of the 
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FIG. 1. — Bj-band image of the Telescopium (AS0851) group with the 
brightest group galaxy NGC 6868 to the left and second brightest group 
galaxy NGC6861 to the right. Other bright members of AS0851 are also 
labeled. Squares show the Chandra ACIS-I fields-of-view (FOVs) for the 
observations discussed in this work. North is up and East is to the left. 



central black hole masses for both NGC 6868 and NGC 6861 
in Table Q] Using the 250 km s" 1 central stellar velocity dis- 
persion for NGC 6868 in the black hole mass - stellar velocity 
dispersion (Mbh-p*) relation, 

10g (^)= (8 - 13± - 06) + (4 - 02 ±0 - 32)10 ^2^oim^)' 

(1) 

where M BH is the mass of the black hole and <r* is the central 
stellar velocity dispersion (Tremaine et al. 2002; Ferrarese 
& Merrit 2000; Gephardt et al. 2000), we infer a central 
black hole mass in NGC 6868 of - 3.3 x 10 8 M Q , typical for 
black hole masses in elliptical galaxies in groups (see, e.g. 
Tremaine et al. 2002, Verdoes Kleijn et al. 2006). An al- 
ternative estimator, the Magorrian relation, relates the mass 
of the central black hole Mbh to the stellar mass M» of the 
galaxy's bulge (as stated in Haring and Rix 2004; see also 
Magorrian et al. 1998; Marconi & Hunt 2003), 

log (^) = (8 - 20±0 - 10)+(1 - 12±ao6)log (lo^)- 

(2) 

For elliptical and spheroidal galaxies the bulge mass is re- 
placed with the total stellar mass of the galaxy. Following 
Gilfanov (2004), we use the 2MASS total K-band luminosity 
of NGC 6868 to estimate the total stellar mass for NGC 6868 
in eq. |2]and find Mbh = 3.1 x 1O 8 M , in excellent agreement 
with the Mbh - &* relation. 

Optical measurements of NGC 6868 reveal a central dust 
lane with weak spiral features (Buson et al. 1993, Veron- 
Cetty & Veron 1988, Hansen et al. 1991), and indicate that 
gas at the center of the galaxy is not moving in regular orbits 
(Zeilinger et al. 1996). IRAS measurements of NGC 6868 at 
60 and 100 /im confirm the presence of cold dust in the galaxy 
(Bregman et al. 1998). A possible explanation for these dust 
features and unusual stellar kinematics is that NGC 6868 has 
recently captured a dust rich companion. Such a cold gas 
rich merger could fuel AGN activity and/or black hole growth 
(Hardcastle et al. 2007; Machacek et al. 2008). NGC 6868 
has been observed in the radio at 2.3, 5, and 8.4 GHz with 
the Parkes and Australian Telescopes (Slee et al. 1994) and 
at 843 MHz in the SUMSS survey (Mauch et al. 2003). The 
galaxy is classified as a low-power, flat spectrum (a ~ 0.07) 
radio source with total flux densities of 139mJy, and 124mJy 
at 843 MHz and 5 GHz, respectively. The brightness tempera- 
ture and spectral slope are inconsistent with thermal emission 



from star forming HII regions, such that the source of the ra- 
dio emission is likely an active galactic nucleus. 

On the other hand, the ~ 400 km s" 1 central velocity disper- 
sion for NGC 6861 , located ~ 206 kpc away from NGC 6868, 
is one of the highest measured for any early type galaxy 
(Wegner et al. 2003; Koprolin & Zeilinger 2000). From the 
Mbh - c* relation (eq. [TJ, we infer a supermassive black hole 
in NGC 6861 of - 2.5 x 1O 9 M , similar to that in M87, the 
central dominant giant elliptical galaxy in the Virgo galaxy 
cluster, and to that in NGC 4649, a giant elliptical galaxy also 
in Virgo (Tremaine et al. 2002). This is an order of magni- 
tude higher than the 2.1 x 1O 8 M black hole mass expected 
for NGC 6861 from the Magorrian relation (eq. 0). 

The dispersion velocity for gas in the inner 30-50pc of 
NGC 6861, close to the black hole sphere of influence, has 
been measured using central emission line widths from Hub- 
ble Space Telescope STIS data (Verdoes Kleijn et al. 2006; 
Beifiori et al. 2009). Both groups confined their analysis to 
galaxies that did not show unusually disturbed stellar kine- 
matics. In each case the measured central gas velocity for 
NGC 6861 is among the highest in their samples, compara- 
ble or exceeding that found in NGC 4486 (M87). However, 
extracting the black hole mass from the measured gas dis- 
persions is highly sensitive to the assumed morphology of 
the gas and, for the thin disk models used by these authors, 
the unknown inclination of the central gas disk. Beifiori 
et al. (2009) argue that the data are consistent with a range 
of black hole masses, i.e. ~ 1.5x 1O 9 M for disk inclination 
angle i = 33° to 3.6 x 1O 8 M for i = 8 1 °, that within the uncer- 
tainties of the scaling relations, span the predictions of either 
the Magorrian or M B h - v* relation. 

These data suggest several possible scenarios. The black 
hole mass for NGC 6861 may be high, as suggested by the 
Mbh~v* relation, and NGC 6861 may itself be the central 
galaxy in a massive subgroup dark matter halo merging with 
the NGC 6868 subgroup. Alternatively, the high central gas 
velocity dispersion in NGC 6861 could be caused, in part, by 
the input of kinetic energy into the gas from the central AGN 
through outflows or turbulence (Verdoes Kleijn et al. 2006; 
Capellari et al. 2000). Then the mass of the central black 
hole needed to explain the high central gas velocities could be 
smaller, falling below the Mbh - prediction and closer to 
that of the Magorrian relation (eq. |2). NGC 6861 does host a 
weak radio source, with core radio power, measured with the 
Australian Telescope, of 6mJy at 5 GHz (Slee et al. 1994). 
The radio source is extended with a core-to-total radio power 
ratio of ~ 0.25. Although there is no evidence for large scale 
radio jets, the 5" (0.7 kpc) angular resolution of the Australian 
Telescope does not rule out substructure on scales of a few 
tens of parsecs nor preclude feedback from previous episodes 
of AGN activity, each of which might heat the central inter- 
stellar medium (ISM), lowering the black hole mass required 
to model the central gas velocity dispersions observed with 
STIS. Finally, the stellar velocity dispersion of NGC 6861 
could be anomalously high due to gravitational interactions 
with NGC 6868 or other group galaxies. In these latter two 
cases, the puzzle is then what interaction history could cause 
such a high stellar velocity dispersion for NGC 6861, while 
not producing obviously unrelaxed stellar kinematics that 
would have excluded the galaxy from the Beifiori et al. (2009) 
and Verdoes Kleijn et al. (2006) samples. 

NGC 6868 and NGC 6861 have been observed in X-rays 
with both the ROSAT and Chandra X-ray Observatories. 
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These observations have been used to measure total X-ray lu- 
minosities and test X-ray scaling laws (Beuing et al. 1999, 
O' Sullivan et al. 2001, Ellis & O' Sullivan 2006) and also 
to determine average density and temperature profiles for the 
galaxies (Fukazawa et al. 2006). In this paper we use the 
Chandra X-ray measurements of the density and temperature 
distribution of hot gas in NGC 6868, NGC 6861, and in the 
AS0851 group core to investigate the interaction history of 
these galaxies and how those interactions may affect nuclear 
activity and black hole growth. Our discussion is organized 
as follows: In §J2] we review the Chandra observations and 
our data reduction and processing procedures. In $3] we dis- 
cuss the mean X-ray surface brightness around each galaxy 
and search for substructure in the AS085 1 group. In $4] we 
analyse density and temperature features observed in these 
data and how they may constrain possible interaction scenar- 
ios between the galaxies and group IGM. In $5] we comment 
on how these interactions may impact the black hole mass re- 
lations. We summarize our results in $6] Unless otherwise 
indicated, quoted uncertainties are 90% confidence levels for 
spectral parameters and la uncertainties for X-ray counts. 
Coordinates are J2000. We adopt a luminosity distance of 
28Mpc, obtained from surface brightness fluctuation mea- 
surements of NGC 6868 (Tonry et al. 2001), as representa- 
tive of the distance to the core of the AS0851 galaxy group 
and to NGC 6861. Assuming the five year WMAP values for 
cosmological parameters (il m = 0.27, tt\ = 0.73, Ho = 0.71; 
Hinshaw et al. 2009), 1" at 28 Mpc corresponds to a distance 
of 0.134kpc in the plane of the sky. 

2. CHANDRA OBSERVATIONS AND DATA REDUCTION 

NGC 6868 was observed for 23.5 ks on 2002 January 11 
(Obsid 3191) and NGC 6861 for 22.8 ks on 2002 July 26 (Ob- 
sid 3190), each with the Advanced CCD Imaging Spectrom- 
eter - Imaging Array (ACIS-I; Gamire et al. 1992, Bautz 
et al. 1998) in VFAINT mode on board the Chandra X-ray 
Observatory. We used the standard X-ray analysis packages 
from CIA03.4 and FTOOLS for image processing and spec- 
tral extraction and XSPEC1 1.3.0 for spectral modeling. The 
data were initially filtered to reject bad grades (1,5,7) and data 
that fell on hot pixels. Data flagged in VFAINT mode as 
having excessive count rates in the border pixels surround- 
ing event islands were also removed to optimize signal to 
noise at energies below 1 keV, where most of the source emis- 
sion is expected for these galaxies. We then reprocessed the 
data and created exposure maps and response files applying 
the most recent gain tables and instrumental corrections, in- 
cluding corrections for the charge transfer inefficiency on the 
ACIS-I CCDs, the time-dependent build-up of contaminants 
on the optical filter, and the secular drift of the average pulse- 
height amplitude for photons of fixed energyQ The data were 
cleaned using a 3er clipping algorithm in the script lc_clean 
to remove periods of anomalously high and low count rates. 
Inspection of the resulting light curves showed that this was 
adequate for OBSID 3191, resulting in a useful exposure of 
21551 s for that observation. However, inspection of the light 
curves for OBSID 3190 (NGC 6861) showed that contamina- 
tion from a strong particle flare (ringing) persisted throughout 
the latter half of the observation. Thus we used only the pre- 
flare data from OBSID 3190 in our analysis, resulting in a 
useful exposure of 9852 s for NGC 6861. 




FIG. 2. — (upper) 0.5 — 2keV Chandra mosaic of diffuse emission in 
NGC 6868 and NGC 6861 in the AS0851 galaxy group. Point sources have 
been excluded and the resulting image has been background-subtracted, ex- 
posure corrected and smoothed with an 8" Gaussian kernel, lpixel = 1 "968 X 
1 ''968. Circles indicate the galaxy centers. North is up and east is to the left. 
(lower) 0.5-2keV image of NGC 6868 and NGC 6861 adaptively binned 
using Voronoi tessellation, with 3"936 X 3"936 pixels and S/N = 5, to high- 
light the fainter structures. 



Backgrounds for the imaging analysis and for spectral anal- 
ysis of the IGM were constructed from the 1 .5 x 10 6 s source 
free background data sefl appropriate for our observation 
dates and instrument configuration. Background data sets 
for each observation were normalized by comparing count 
rates in the 9.0- 11.5keV energy band, where particle back- 
ground dominates. This resulted in additional normalization 
factors of 0.98 and 1 .07 for the backgrounds for OBSID 3191 
and OBSID 3190, respectively. We identified X-ray point 
sources in each observation in the 0.3-7keV energy band 
using a wavelet decomposition algorithm with a 5a detec- 
tion threshold. Using 10 photon counts, we identify point 
sources to limiting 0.3-7 keV fluxes of 8.6 x 10~ 15 erg s" 1 and 
4.3 x 10" 15 ergs" 1 for OSID 3190 (NGC6861) and OSID 
3191 (NGC 6868), respectively. Since we are interested in the 
properties of the diffuse X-ray emission and possible AGN 
activity, X-ray point sources, other than the galaxies' nuclei, 
were removed from all subsequent analyses. 

3. MEAN SURFACE BRIGHTNESS PROFILES AND 
GROUP STRUCTURE 

In the upper panel of Figure [2] we show a mosaiced 0.5- 
2keV Chandra image of the diffuse emission surrounding 
NGC 6868 (left) and NGC 6861 (right) in the galaxy group 
AS085 1 . The image has been background subtracted and ex- 
posure corrected, and, after exclusion of X-ray point sources, 
has been smoothed with an 8" Gaussian kernel. In the lower 
panel of Figure [2] we present an adaptively binned 0.5-2keV 
Chandra image of the same regions, using the Voronoi tessel- 
lation technique with a signal-to-noise {S/N) of 5 to high- 
light faint emission features. We first construct mean sur- 
face brightness profiles for NGC 6868, NGC 6861 and the 
surrounding gas to understand the structure of the group as 



see, http://cxc.harvard.edu/contrib/alexey/tgain/tgain.html 



2 see http://cxc.harvard.edu/contrib/maxim/acisbg 
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FIG. 3. — {upper) Azimuthally averaged 0.5 — 2keV radial surface bright- 
ness profile of NGC 6868 constructed using concentric circular annuli, cen- 
tered on NGC6868's nucleus, with logarithmically increasing width. The 
solid line denotes the sum of the double /3-model shown in the lower panel 
and a constant (7.5 X 10~'° photons cirr 2 arcsec~ 2 , dotted line) residual soft 
Galactic X-ray background, (lower) The azimuthally averaged surface bright- 
ness profile of NGC 6868 from the upper panel after subtraction of the soft 
Galactic X-ray component (dotted line). The solid line is the sum of two 
/3-models (shown as dashed lines), an inner model with index = 0.65 and 
core radius r c = 0.492 kpc and an outer model with /3 = 0.52 and r c = 1 1 kpc. 
The dotted line denotes the constant soft background level and thus when 
signal-to-noise S/N = 1. 



FIG. 4. — (upper) 0.5 — 2keV Chandra X-ray image of the bright inner 
emission region of NGC 6861 with 2MASS K-band contours superposed. Re- 
solved point sources have been excluded and the image has been smoothed 
with a a = 0"984 Gaussian kernel. K band contour levels are 560, 580, 610, 
720 data units, respectively. The dashed arrow shows the 2MASS K band 
position angle (-40° counter-clockwise from North), while the solid arrow 
is the position angle (—52° ±7°) determined from the elliptical /3-model fit 
(/3 = 0.62,r c = 0.62 kpc, e = 0.4) to the 0.5 - 2 keV X-ray surface brightness 
profile within 7.3 kpc (54") of NGC 6861. lpixel = 0"492 x Of'492. (lower) 
Azimuthally averaged surface brightness profile of NGC 6861 to large radii. 
Solid line denotes the double spherical /3-model fit where /3 = 0.61(0.38) and 
r c = 0.5(1 1 .9) kpc for the inner (outer) /3-model component, respectively.) 



a whole and search for evidence of substructure in AS0851 
that might signal the merger of two subgroups. In the follow- 
ing section we will discuss the evidence for interactions 
suggested by the surface brightness features identified in Fig- 
ure |2]in and around each galaxy. 

3.1. NGC 6868 

In Figure [3] (open circles, upper panel) we show the mean 
radial surface brightness profile for NGC 6868 constructed us- 
ing concentric circular annuli of logarithmically increasing 
width, centered on NGC6868's nucleus. The profile flattens 
at large radii indicative of a residual soft X-ray background 
component in the data. ROSAT All Sky Image (RASS) maps 
confirm that NGC 6868 and NGC 6861 lie in a region of en- 
hanced soft Galactic X-ray emission. Since this soft Galac- 
tic X-ray background is expected to be constant over the 
field of view, we estimate this background using the mean 
0.5-2keV emission in the S2 CCD from OBSID 3190, 
the most distant (r mean ~ 180 kpc) source free region from 



NGC 6868, where we expect the soft Galactic background 
to dominate. We find the soft Galactic background level to 
be 7.5 x 10~ 10 photons cm" 2 arcsec" 2 and denote this level as 
the horizontal dotted line in both panels of Figure [3] In the 
lower panel of Figure [3] we show the surface brightness pro- 
file for NGC 6868 in the 0.5-2keV energy after subtraction 
of this residual soft Galactic background. The mean sur- 
face brightness profile is well parameterized by the sum of 
two /3-models, an inner /3-model with /3 = 0.65 and core ra- 
dius r c = 0.49 kpc, and an outer /3-model with j3 = 0.52 and 
r c = 11 kpc. Note that beyond r ~ 35 kpc the fit becomes 
highly uncertain since the soft Galactic background domi- 
nates the faint group emission in this energy band. 

3.2. NGC6861 

The distribution of starlight in NGC 6861 is highly eccen- 
tric (see Fig. Q]). From 2MASS K band measurements, the 
ellipticity (e = 1-b/a) within the "total" K-band isopho- 
tal semi-major radius a = 97 ."4 (13 kpc) is 0.4 with posi- 
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tion angle, measured counter-clockwise from north (J2000) 
of 140° (NED). The brightest X-ray emission, shown in Fig- 
ure!?] is a ls° eccentric and highly compact. Using the CIAO 
tool Sherpa to fit a single two-dimensional elliptical /3-model 
plus constant background to this central emission for r < 54" 
(7.3 kpc), we find an eccentricity 0.40 ± 0.04, consistent with 
that of the starlight, but with position angle 128° ±7°, twisted 
by less than 2a relative to position angle measured for the 
stellar light. We find j3 = 0.62^° jjj and a core radius r c = 
0.62^03 kpc, similar to elliptical beta model fits to other el- 
liptical galaxies in galaxy groups (Osmond & Ponman 2004). 
As we increase the fit radius to include fainter emission, the 
eccentricity decreases (<? = 0. 19 ± 0.05 for r < 42 kpc), as ex- 
pected as we include emission from gas residing in the more 
spherical gravitational potential of the group or subgroup. To 
separate the possible subgroup emission from the galaxy com- 
ponent, we fit the background-subtracted, azimuthally- aver- 
aged surface brightness profile, shown in Figure |4] with the 
sum of two sp heric ally symmetric beta models (as we did for 
NGC6868 in f3~TT >. The soft Galactic background denoted by 
the dashed line in Figure|4]has also been subtracted. We find 
(13, r c ) of (0.61, 0.5 kpc) for the inner (galactic) component and 
(0.38, 11. 9 kpc) for the extended (associated subgroup) com- 
ponent. The extended component j3 parameter for NGC 6861 
is much shallower than that found for NGC 6868, but still falls 
within the range of j3 parameters (0.36 < j3 < 0.58) found for 
the extended X-ray components around elliptical galaxies in 
groups in the GEMS survey (Osmond & Ponman 2004). 

3.3. Group Structure 

The fact that the surface brightness distribution in and 
around NGC 6868 and NGC 6861, the two brightest galax- 
ies in AS0851, are both described by double /3-models sug- 
gests that these galaxies may be the dominant galaxies in 
two subgroups that are merging. To investigate whether there 
is evidence for subgroup substructure in the radial veloc- 
ity distribution of galaxies in the AS0851 group, we iden- 
tify all galaxies with measured radial velocities 2000 < v r < 
4000 km s" 1 within 4° (2Mpc) of the brightest group galaxy 
NGC 6868 (NED; Huchra & Geller, CfA Redshift Survey 
2008)0 We find 46 galaxies, with radial velocities closely 
grouped between 2400 < v r < 3700 km s" 1 and well removed 
in redshift from prominent, background large scale filaments 
at v r ~ 5000 km s" 1 that appear nearby in projection (see the 
upper left panel of Fig. |5). Using the method of Osmond 
& Ponman (2004), we find a velocity dispersion for these 46 
galaxies of 274 ± 29 km s" 1 . While this velocity dispersion 
is typical of that expected for a cool 1 keV group, the spa- 
tial extent is much larger, i.e. is expected to be only 
~ 0.6 Mpc for a IkeV group (Osmond & Ponman 2004), 
which may argue in favor of the merger of two cool galaxy 
subgroups. However, as seen in the lower left panel of Fig- 
ure [5] there is little evidence for substructure in radial ve- 
locity space even on these larger scales. 65% of the galax- 
ies in this sample (30 of 46) are clustered within 1 Mpc of 
NGC 6868 or NGC 6861 and 52% (24 galaxies) are within 
rsoo (0.6 Mpc) of NGC 6868. The group velocity dispersion of 
galaxies within 2 Mpc, 1 Mpc and 0.6 Mpc (274 ± 29 km s" 1 , 
263 ± 34kms _I , and 251 ± 37kms, respectively) all agree 
within errors. NGC 6868 and NGC 6861 have radial veloci- 
ties separated by only 35kms _1 and both radial velocities lie 

3 http://www.cfa.harvard.edu/~huchra/zcat/. 



close to the median radial velocity of the galaxy sample as 
a whole (~ 2853 km s _1 ). Thus the merger of the NGC 6868 
and NGC 6861 subgroups is predominantly in the plane of the 
sky, and it is unclear which subgroup, if either, may be at the 
center of the AS085 1 group potential. 

4. EVIDENCE FOR INTERACTIONS 

Galaxy interactions and the dynamical evolution of the sur- 
rounding IGM are revealed through asymmetries in the den- 
sity and temperature of the gas. When the relative motion 
between the galaxy or subgroup gas halo and any ambient 
medium lies close to the plane of the sky, we expect ram 
pressure acting on the infalling body to produce a sharp X- 
ray surface brightness discontinuity along the leading 'edge' 
of the infall trajectory, associated with a cold front, as well 
as a trailing tail or wake of stripped galaxy or subgroup gas 
(see, e.g. Forman etal. 1979, Rangarajan etal. 1995, Randall 
etal. 2008 forM86; Scharf et al. 2004, Machacek et al. 2005a 
for NGC 1404; Machacek et al. 2006 for NGC 4552). The 
surface brightness morphology of the gas around NGC 6868 
suggests motion. In particular, the surface brightness distribu- 
tion around NGC 6868, shown in Figure [2] appears flattened 
to the north, with a sharp discontinuity (edge) to the north 
and west. X-ray emission from NGC 6868 is more extended 
(tail-like) to the south and east. Other asymmetries in the sur- 
face brightness distribution include a bright clump of emis- 
sion ~ 1' (8 kpc) to the east of the nucleus, and a roughly 
spherical deficit in surface brightness in the IGM to the north- 
west of NGC 6868, as well as a partial rim-like feature in the 
ISM gas to the southwest. These latter features may be evi- 
dence for X-ray cavities inflated by prior AGN activity. 

Although the brightest X-ray emission in NGC 6861 fol- 
lows the optical isophotes, it is extremely compact for a mas- 
sive galaxy and the asymmetries in the X-ray surface bright- 
ness revealed at larger radii are striking (see Figs. [2] and 
[6]). An approximately 3'(24kpc) long X-ray tail (W Tail) ex- 
tends due west of NGC 6861. A plume of emission rises 70" 
(9.4) kpc to the northeast, perpendicular to the major axis of 
the galaxy, and a tail of X-ray emission trails NGC 6861 to 
the northwest (NW tail). 

Several scenarios offer possible explanations for the dis- 
torted surface brightness morphologies observed in and 
around NGC 6868 and NGC 6861. First, NGC 6868, 
NGC 6861 or both galaxies may be moving with respect to 
the group IGM, such that the extended tail-like features may 
be the result of ram pressure (Gunn & Gott 1972) and/or 
turbulent-viscous (Nulsen 1982; Quilis et al. 2000) stripping 
of galaxy gas, or gravitational focusing of the group IGM 
(Bondi & Hoyle 1944; Bondi 1952; Ruderman & Spiegel 
1971; Machacek et al. 2005b). Second, we may be wit- 
nessing a more evolved merger of two massive subgroups 
whose dominant galaxies are NGC 6868 and NGC 6861, re- 
spectively, such that much of the NGC6861's outer stellar 
halo and subgroup gas has been tidally stripped after multiple 
passes through the core of AS0851 (Byrd & Valtonen 1990; 
Moore et al. 1996) and the group IGM set in motion and 
mixed as a result of these gravitational encounters (sloshing, 
see, e.g., Markevitch et al. 2001; Mazzotta etal. 2003; Tittley 
& Henriksen 2005; Ascasibar & Markevitch 2006). Third, the 
disturbed morphologies may partly be the result of energetic 
outbursts from the active nuclei at the centers of the dominant 
group galaxies, possibly triggered by merger activity (see, e.g. 
David et al. 2009 for a study of the NGC 5044 group). The 
dynamical processes in these scenarios may act individually 
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FIG. 5. — (upper left) Large scale structure environment of the Telescopium (AS0851) group from the CfA Redshift Survery 2008. Colors denote radial velocity 
bins of < 1500 km s" 1 (magenta), 1500 -2999 km s" 1 (red), 3000- 4499 kms(blue), 4500 -5999 km s _1 (green). The arrow denotes the position of the AS0851 
galaxy group. Shapes indicate galaxy morphology based on T type: solid circles (ellipticals, T = -7 to— 3); open circle (spheroidals, T = -2 to 0); asterisk ( early 
type spirals Sa-c, T = 1 to 5); crosses (late type spirals and irregulars Sd - Im, T = 6 to 10); triangles (unclassified) and increasing symbol size denotes decreasing 
magnitude. Coordinates are J2000.0. (upper right) Galaxy distribution within 2° of the brightest group galaxy NGC 6868 with the same symbols as in the upper 
left panel, (lower left) Radial velocity histogram of all galaxies with radial velocities 2000 < v,- < 4000 km s within 4° (2Mpc) of NGC 6868. (lower right) 
Radial velocity histogram of all galaxies with radial velocities 2000 < v r < 4000 km s within a 1 Mpc (124') radius of NGC 6868. 



or together to effect the transformation of these galaxies and 
their environment. 

4. 1 . NGC 6868 's Edge and Tail: the Case for Infall 
4.1.1. Modeling the Edge 

The sharp X-ray surface brightness edge to the north of 
NGC 6868 and trailing tail to the south (see Fig. [2), sug- 
gest that NGC 6868 may be moving with respect to the group 
IGM. In Figures [7J and [8] we investigate the northern edge- 
like feature. Following Vikhlinin et al. (2001), we construct 
the 0.5-2keV surface brightness profile for NGC 6868 across 
the edge using elliptical annuli, centered on NGC 6868 's nu- 
cleus and constrained to lie in an angular sector extending 
from 41.5° to 164.5°, to avoid the northwestern cavity. The 
elliptical annuli are chosen concentric to a bounding ellipse 
that traces the surface brightness discontinuity in the sector 
of interest (see Fig. [7ji. The width of each elliptical annulus 



increases (decreases) with constant logarithmic step size as 
one moves outward (inward) from the bounding ellipse. Be- 
yond 50 kpc, the elliptical annuli are no longer fully contained 
within the ACIS-I field-of-view. We define three narrow rect- 
angular regions (also shown in Fig. [7} to extend the surface 
brightness profile to larger radii. 

The resulting profile is shown in Figure [8] In the upper 
panel of Figure [8] we compare the surface brightness profile 
across the northern edge (filled squares) to the azimuthally 
averaged profile for the galaxy (open circles). We see a sharp, 
statistically significant, drop in surface brightness between 
20 < r < 35 kpc, consistent with the presence of an edge. To 
model the X-ray surface brightness profile across the edge, we 
integrate the square of the electron density n e times the corre- 
sponding X-ray emissivity A along the line of sight, assuming 
that the electron density on either side of the discontinuity is 
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FIG. 6. — (upper) Background subtracted, exposure corrected 0.5-2keV 
Chandra X-ray image of NGC6861 showing a plume of emission to north- 
east, a bright tail to the west, and a fainter tail-like feature extending to the 
northwest. The image has been smoothed with an 8" gaussian kernel, (lower) 
Projection of surface brightness in the rectangular region in the upper panel, 
showing a factor ~ 2 enhancement in surface brightness in the northwestern 
tail over that in the adjacent IGM. 



well described by spherical power law models of the form: 



n c = «i(— — ) 



f edg 
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where r ec jge is the location of the edge. The observed surface 
brightness discontinuity is then given by 
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where A,- and n -, are the X-ray emissivity and electron density 
for gas inside (z = 1) and outside (z = 2) the edge. We fit the 
surface brightness profile across the edge using a multi-variate 
chi-square minimization scheme allowing the edge position 
'"edge, electron density power law index ot\ inside the edge, 
and the size of the discontinuity (jump) J to vary. Since the 
signal-to-noise drops below 1 just outside the edge (see Fig. 
[8]), the powerlaw index a2 for the IGM gas density is not well 
determined by these data. We fix ct2 = -1.56 from the pa- 
rameterization of the mean surface profile at large radii. The 
results of our model fits are listed in Tableland plotted in the 
lower panel of Figure [8] where we have varied the fit range 




FIG. 7. — 0.5-2keV blank-sky background subtracted, exposure corrected 
Chandra image of NGC 6868 with the angular sector (41.5° to 164.5°), 
bounding ellipse centered on NGC6868's nucleus with (major, minor) axes 
and position angle of (228" 3, 131"9) and 96°, respectively. Outer box re- 
gions, used to construct the surface brightness profile at large radii, are su- 
perposed. 



to probe the sensitivity of the edge position and jump discon- 
tinuity to the low S/N data. We find r e( [ge ~ 23 ± 3 kpc and 
jump discontinuity J ~ 2.2^ j. 

To complete the density model, we compute the X-ray 
emissivity A in the 0.5-2keV energy band from APEC ther- 
mal plasma models for the X-ray emission from elliptical an- 
nular sectors, listed in Table [3j chosen concentric to the edge 
bounding ellipse (see Fig. |7]l inside (NEdge_i) and outside 
(NEdge+o,NEdge+i) the edge. The results of these spectral fits 
are given in Table [4] Although the uncertainties are large, the 
temperature of the IGM outside the edge to the north is higher 
(0.73^;g keV) than inside the edge (~ 0.64 ± 0.08 keV), sug- 
gesting that the edge is a cold front. For temperatures < 1 keV, 
the Fe L complex line emission dominates the thermal con- 
tinuum in the 0.5-2keV energy band, such that the X-ray 
emissivity depends sensitively on the abundance of the gas, 
and can vary by a factor ~ 5 as the abundance varies from 
0.1 to 1Z Q . This introduces an uncertainty in the inferred 
density ratio across the edge by a factor ~ 2. Although the 
abundances are poorly constrained by these data, the spectral 
models in Table |4] suggest that low abundances (A < 0.3 Z©) 
are favored in the group IGM north of the edge, consistent 
with abundances measured in other cool groups and in the out- 
skirts of clusters (Rasmussen & Ponman 2007; Finoquenov 
et al. 2000; De Grandi et al. 2004). The spectral models in 
Table |4] also suggest that higher abundances A ~ 0.5- 1Z© 
are favored inside the edge closer to the center of NGC 6868, 
consistent with abundances measured in other dominant ellip- 
tical galaxies (see, e.g. Finoguenov & Jones 2001; Matsushita 
et al. 2003). 

For simplicity we define a pair of fiducial spectral models 
with temperature kT = 0.73 keV and abundance A = 0.3Z© 
for IGM gas to the north, outside the edge, and kT = 0.64 keV 
and A = 0.5Z© for galaxy gas inside the edge. Normalizing 
the density model using the IGM gas density from eq. [3] at 
r = 34 kpc to the north, where the contribution from the galaxy 
gas is negligible (see Fig. [3]), we find an electron density in 
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FIG. 8. — (upper) 0.5 — 2keV surface brightness profile across the northern 
edge of NGC6868 (filled squares), using the regions shown in[7j superposed 
on the azimuthally averaged 0.5-2keV surface brightness profile from the 
upper panel of[J](open circles). The lines are defined as in Fig. [J] (lower) 
Model fits to the northern surface brightness edge in NGC 6868 in the sector 
shown in Fig. [7J Open squares denote the surface brightness profile after 
subtraction of the soft Galactic background. Curved lines denote the density 
model fits to the background- subtracted surface brightness profile across the 
edge over the radius range 7.9 -34 kpc (solid black), 7.9 — 68 kpc (dot-dash 
blue), 5.8 — 35 kpc (long dash violet), and 5.8 — 68 kpc (short dash red) listed 
in Tableff] The corresponding vertical lines denote the edge position in each 
case and the dotted horizontal line shows where S/N = 1. 



the IGM n e (r = 34 kpc) = 4.4 x 10" 4 cm" 3 . Combining this 
with the density model fits to the edge location and surface 
brightness jump from eq. |4j we find electron gas densities 
at r = r ec [ge inside and outside (ni) NGC6868's north- 
ern surface brightness edge to be (1.2— 1.5) x 10~ 3 cm~ 3 and 
(0.66-1.0) x 10" 3 cm" 3 , respectively (see Table©. Note that 
if there is no abundance gradient across the edge, the density 
ratio derived from eq. [4] would be 23% larger. 

4.1.2. Velocity Constraints 

To constrain the relative motion of NGC 6868 with respect 
to the ambient IGM, we assume uniform gas flow past a rigid 
spheroid, such that the difference between the gas pressure in 
the undisturbed IGM (free stream region) and that at the stag- 
nation point of the flow, where the relative velocity of the IGM 
gas and galaxy gas vanishes, is the ram pressure due to the 
motion of the galaxy through the IGM (Vikhlinin et al. 2001). 
The pressure ratio po / pi between IGM gas at the stagnation 
point and that in the free-stream region allows us to deter- 



mine the Mach number of the flow (Landau & Lifshitz 1959) 
and thus the three-dimensional velocity of the infalling galaxy 
with respect to the ambient gas. Since the stagnation region 
is small and difficult to observe, and pressure is continuous 
across a cold front edge, we use the pressure in NGC 6868 
just inside the northern edge (pi) as a proxy for the pressure 
at the stagnation point. We use the density and spectral model 
fits to gas inside and outside NGC 6868 's northern edge, given 
in Tables [2] and [4] to compute the range of pressure ratios 
(pijpi) and Mach numbers consistent with the fits to the sur- 
face brightness discontinuity (see Fig. [8j. Our results are 
given in Table [5] 

For all models, the inferred velocity for NGC 6868 rela- 
tive to the surrounding group gas is at most transonic, i.e. 
v < 475 kms (Mach < 1.08), given the speed of sound c s = 
440^30 km s" 1 for 0.7 keV group gas. This is similar to that 
found in other infalling galaxies, such as NGC 1404 in For- 
nax (Machacek et al. 2005a) and, if we assume NGC 6861 
rather than NGC 6868 is at the group center, not unexpected 
given the ~ 200 kpc (~ 0.3r5oo) separation between the two 
galaxies. On the other hand, the uncertainties in the ob- 
served metal abundances are large, and, since the spectrum 
for cool (< 1 keV) gas is dominated by the Fe L complex 
line emission, the emissivity of X-ray gas is highly depen- 
dent on the metal abundance of the gas. Solar abundances are 
often measured in the centers of dominant elliptical galaxies 
in groups (for a compilation, see, e.g. Rasmussen & Ponman 
2007 and references therein). If we assume that solar abun- 
dances extend all the way to the surface brightness disconti- 
nuity in NGC 6868, the sharp jump in surface brightness may 
be explained by the ratio of abundances (and thus emissivi- 
ties) across the edge alone. In that case, NGC 6868 may be 
at rest with respect to the surrounding IGM and likely at the 
center of the group potential. However, then stripping could 
not explain the origin of the southern tail. 

4.1.3. NGC 6868' s Southern Tail 

As suggested in Figure [2] and confirmed in the lower panel 
of Figure [TT] X-ray emission (and thus gas density) is higher 
to the south of the galaxy (171.5° < <fi < 351.5°) than to the 
north at the same radii. If the southern tail is composed of 
ram-pressure stripped galaxy gas, we would expect the metal 
abundance in the southern tail to be higher, similar to abun- 
dances inside the edge in the galaxy gas halo, than in the IGM 
to the north. To test whether the southern tail is gas ram- 
pressured stripped from NGC 6868, we fit the spectrum from 
an elliptical annular sector (S+\) at the same mean radial dis- 
tance and with the same radial bin width as sector NEdge+i 
to the north (see Table [3), but constrained to lie in the tail to 
the south between 213° to 327.5°. As shown in Table|U the 
best fit single APEC model for this southern region gives a gas 
temperature and abundance of 0.63^og keV and O.29^o^gZ , 
respectively. While the temperature of the gas in this region is 
lower than that of the IGM to the north and comparable to the 
temperature of galaxy gas within NGC 6868, as one might ex- 
pect if the gas were ram-pressure stripped galaxy gas, the un- 
certainties in the measured abundance are too large to defini- 
tively determine the origin of gas in the southern tail. Deeper 
X-ray data are needed to distinguish between: (1) cool, higher 
abundance (~ 0.5 Zs) ram-pressure stripped galaxy gas, (2) 
lower temperature, lower entropy subgroup gas with an abun- 
dance comparable to that found in the IGM to the north, or (3) 
a multiphase mixture of the two, possibly signaled by the low 
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FIG. 9. — Gas Temperature as a function of mean photon energy in the 0.7 — 
1.2 keV enery band for ACIS-I CCD 10 (circles), II (squares), I2(triangles), 
and 13 (astericks), assuming a single temperature APEC model with fixed 
metal abundance A = 0.5Zq and Galactic absorption of 3.9 X 10 20 cm -2 . 
The solid line denotes the best linear fit (kT = -6.258 + 7.5492£) restricted 
to0.4<itr< 1.3 keV. 



90% CL lower bound on the abundance 0.11 Z©) from the 
single temperature APEC model fit (Buote 2000). 

4.2. Temperature Asymmetries: the Case for Sloshing 

Modest density jumps and subsonic relative velocities, as 
observed for NGC 6868, are also characteristic of gas 'slosh- 
ing', where the surrounding IGM has been set in motion by 
the gravitational perturbation of a prior passage of galaxies or 
galaxy subgroups past each other in the group core (Marke- 
vitch etal. 2001). Simulations of sloshing in subcluster merg- 
ers show that the signatures of sloshing are most pronounced 
in the temperature structure of the gas, often producing cool, 
spiral-shaped tail-like features, as higher entropy gas from 
larger radii is mixed with the nearby lower entropy galaxy 
and subgroup gas (Ascasibar & Markevitch 2006). 

4.2.1. Temperature Map 

Since NGC 6868 and NGC 6861 are elliptical galaxies in 
only a moderately massive galaxy group, we expect tempera- 
tures for gas in the galaxies and surrounding group/subgroup 
gas to be < 1 keV. For such low gas temperatures the X-ray 
spectrum peaks at ~ 1 keV and is dominated by lines in the 
Fe L complex. Using XSpec simulations of single temper- 
ature APEC plasma models, we find that there is nearly a 
linear correspondence between the temperature of the X-ray 
emitting gas in the 0.4- 1.3keV range and the mean photon 
energy measured in the 0.7- 1.2keV energy band, contain- 
ing the Fe L emission (See Fig. [9] and David et al. 2009). 
We exploit this correspondence to construct the tempera- 
ture map of NGC 6868 and NGC 6861, shown in Figure [10] 
After resolved point sources have been removed, we adap- 
tively smooth the narrow (0.7- 1.2 keV) band Chandra sur- 
face brightness images for NGC 6868 and NGC 6861 to deter- 
mine regions of interest. We then average the photon energy 
in the Fe L peak over these regions and use the calibration 
shown in Figure [9] to relate these mean energies to the gas 
temperature. To check that our results were not an artifact of 
the binning procedure, we constructed a second temperature 
map, shown in the lower left panel of Figure [10] by explicitly 
fitting single temperature APEC models to regions containing 



a minimum of 500 counts (O' Sullivan et al. 2005; Maughan 
et al. 2006; Randall et al. 2008). Although the second bin- 
ning is much coarser and averages over the nuclear region in 
each galaxy, the results for the extended emission outside the 
nucleus are qualitatively the same. 

The temperature maps confirm that the gas temperature in 
NGC 6868 inside (south of) the surface brightness edge is 
lower than in the group IGM to the north outside the edge. 
The lowest temperature gas (~ 0.59 keV) in NGC 6868 is dis- 
placed - 10-20 kpc (i:2-2.'5) to the east of NGC6868's 
nucleus. There is a ring of higher temperature gas close 
to the nucleus, possibly indicating nuclear activity. The ex- 
tended gas tail to the south of NGC 6868 is cool (~ 0.65 keV), 
and appears to spiral to the south and southwest around a 
~ 30 kpc long finger of higher temperature emission south- 
west of the nucleus. Both temperature maps show the spi- 
ral feature in NGC6868's extended southern tail. Although 
the temperature of the southern tail is consistent with that 
expected for ram-pressure stripped galaxy gas and such tails 
may curve, following the orbit of the infalling galaxy or sub- 
group, the tightly wound spiral pattern of cool gas wrapping 
about a higher temperature filament is more characteristic of 
the 'sloshing' signatures seen in simulations of IGM gas set in 
motion by interactions of galaxies or subcluster cores near the 
centers of clusters (Tittley & Henriksen 2005; Ascasibar & 
Markevitch 2006). Such features have been observed recently 
in X-ray observations of sub-cluster mergers and in galaxy 
groups (e.g see Markevitch & Vikhlinin 2007 for A2142; 
Mazzotta & Giacintucci 2008 for MS 1455.0+2232; Johnson 
et al. 2009 for A1644; Randall et al. 2009 for the NGC 5098 
group). 

The temperature structure of NGC 6861 is dramatically dif- 
ferent. In both temperature maps, single temperature APEC 
models suggest that the temperature of gas in the bright west- 
ern tail is higher than in the surrounding IGM. Such high tem- 
perature tails have been observed in other galaxy groups, but 
are rare (see, e.g. Osmond et al. 2004 for the NGC5171 
group; Machacek et al. 2005b for the Pavo group). The upper 
panels of Figure [10] also show a sheath of high temperature 
gas ~ 8 kpc (1') from NGC6861's center, well outside the X- 
ray bright central region of NGC 6861. There is no feature in 
the temperature maps corresponding to the northwestern tail 
observed in the X-ray surface brightness images (Figure©. 

4.2.2. Spectral Modeling 

We extracted spectra from regions based on the temperature 
map features for detailed fitting using XSpec 1 1 .3.0. These re- 
gions are overlaid on the temperature and surface bri ghtn ess 
maps in the upper right and lower right panels of FigureTlOl re- 
spectively, and are listed in Table|6]for NGC 6868 and Table[8] 
for NGC 6861. The data were binned to ensure a minimum of 
20 counts per bin and also using constant logarithmic bin size. 
We found no significance differences in the results obtained 
from these two binning methods. Backgrounds were taken 
from rescaled blank sky background sets (see §|2), extracted 
from the same region as each source. We model the spec- 
tra using absorbed APEC thermal plasma (Smith et al. 2001 
) or APEC plus power law models and fit the data over the 
0.5-7keV energy range. To check our sensitivity to the soft 
Galactic background, we refit the data over the 0.7-7keV 
range where this soft Galactic background will not contribute 
and found the fits unchanged. 
NGC 6868: 

For NGC 6868 we concentrate on features in and to the 
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FIG. 10. — (upper left) Temperature map of NGC 6868 and NGC 6861 constructed from the mean photon energy in the Fe peak at 0.7 - 1.2 keV. Rectangular 
regions correspond to fields fit with the coarse temperature map shown below in the lower left panel (upper right) Same temperature map with spectral regions 
from Tables [6] and [8] superposed. N and G mark the centers of NGC 6868 and NGC 6861, respectively, (lower left) Coarse temperature map constructed by 
explicitly fitting single temperature APEC models to regions containing a minimum of 500 counts, (lower right) Spectral regions from Tables [6] and [8] on the 
adaptively binned 0.5 - 2 keV Chandra X-ray image of NGC 6868 and NGC 6861 from Fig. |2] 



south of the galaxy. We extract spectra for the nuclear re- 
gion of the galaxy (N), the lowest temperature region east 
of NGC6868's center (cold clump C), the filament of higher 
temperature gas to the southwest (F) and the lower tempera- 
ture, spiral tail-like feature (the combined regions T). The re- 
sults of our spectral fits for NGC 6868, listed in Table|7] gen- 
erally confirm the results of Figure Q7J The galaxy nucleus 
is well described by a two component APEC + power law 
model with gas temperature 0.63^° "gkeV, solar abundance, 
and power law index ~ 1 .4 to account for unresolved X-ray 
binaries in the central region of the elliptical galaxy and/or 
a weak AGN. Outside NGC6868's central region, the spec- 
trum of the higher temperature filament (region F) and the 
southern spiral-like tail (region T) are both well described 
by single APEC thermal plasma models with temperatures of 
0.70 ±0.05 and 0.62 ±0.04 keV, respectively. APEC models 
for the filament F have a modestly improved reduced \ 2 for 
low abundances (~ 0.3-0.2Z Q ), suggesting that the filament 
may be composed of IGM rather than galaxy gas. The most 
problematical region to interpret in the temperature map of 
NGC 6868 is the cool clump region C. While a statistically 
acceptable fit (null = 0.09) is found for a single temperature 
(kT = 0.59+f] :E)5 keV ) APEC model, there are significant resid- 
uals in the 0.7- 1 keV region and the x 2 /dof = 38/28 is not 
particularly good. Addition of a second APEC or power law 
component did not improve the fit. 
NGC 6861: 

To determine the thermodynamics properties of gas sur- 
rounding the spheroidal galaxy NGC 6861, we measure the 
gas temperature in three rectangular regions: region E lo- 
cated 43.2 kpc to the east of NGC6861's nucleus in the di- 
rection of NGC 6868, region NW located to the northwest 
coincident with the northwestern tail shown in Figure [6] and 



in the lower right panel of Fig. QJJ and region WT contain- 
ing the bright tail to the west of NGC 6861. We also model 
spectra within NGC 6861 using an elliptical region G centered 
on NGC6861's nucleus with semi-major (minor) axes of 4.2 
(2.6) kpc, respectively, containing the nucleus and brightest 
galactic emission (see Fig. HJ, and in a forward annular sector 
FA, with inner and outer radii of 5.3 kpc and 14.3 kpc, respec- 
tively, and constrained to lie between 80° and 282° measured 
counter-clockwise from west, containing the plume (Fig. [6) 
and 'sheath', identified in Figure[l0] surrounding the galaxy's 
central region. The results of our spectral fits to these regions 
are listed in Table [9] The central region of the galaxy (G) is 
well modeled by a 0.66^ ;^ keV APEC model, describing the 
contribution of the diffuse galaxy gas, plus a power law com- 
ponent with photon index 1.86 ±0.6, consistent with an ac- 
tive nucleus or population of unresolved X-ray binaries. The 
temperature of gas in the forward annulus region (FA) is sig- 
nificantly higher (1.27 ±0.15keV), in agreement with Fig- 
ure [10] and higher metal abundances (> 0.5Z Q ) are favored, 
as expected for gas associated with the galaxy. Gas in the 
western tail is well described by a single temperature APEC 
model, with gas temperature kT = l.lS^JSj and abundance 
O.25+0J3Z©, confirming that the temperature in the western 
tail (WT) is also high. The low best fit metallicity in region 
WT, albeit with large uncertainties, favors a large IGM com- 
ponent rather than a purely galaxy origin for the gas in the 
western tail. 

The properties of the X-ray emission to the northwest and 
east of NGC 6861 are more uncertain. The temperature and 
abundance of gas to the northwest of NGC 6861, along the 
possible northwestern tail, are 0.90!j];f keV and O.22^;?fZ , 
with similar abundance but lower temperature than gas in the 
western tail. In contrast, a single temperature APEC model 
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fit to the spectrum in region E east of NGC 6861 (denoted fit 
El in Tabled gives a low gas temperature (~ 0.62^ [^keV), 
consistent at the 90% CL with the gas temperature to the north 
of NGC 6868. However, the best fit abundance is anoma- 
lously low (< 0.06Z©) in that model when compared to metal- 
licities outside dominant galaxies in other galaxy groups, 
and the residuals show a 'see-saw' pattern often associated 
with multi-temperature gas (Buote 2000). A two-component 
APEC model (denoted E2 in Table [9]) with a more realistic 
abundance value of 0.2Z Q , consistent with the best fit abun- 
dances found in the northwest and western tail regions, yields 
gas temperatures of ~ 0A6^\lkcV and ~ 1 .24+jj ^ keV, but 
the reduced \ 2 is unchanged over that for the single APEC fit. 
Deeper observations are required to disentangle the complex 
temperature structure to the east of NGC 6861. 

4.2.3. NGC 6861's Tails: Stripping or Sloshing? 

Since bifurcated tails are common features in ram-pressure 
stripped galaxies, such as those infalling into the Virgo clus- 
ter (see, e.g. Randall et al. 2008), NGC 686 1 's bright western 
and fainter northwestern tails might suggest that it is infalling 
towards NGC 6868 near the core of the AS085 1 galaxy group. 
However, most ram-pressure stripped tails are composed of 
galaxy gas with higher metallicities and lower gas tempera- 
tures than the surrounding group IGM, in contradiction with 
the > 1 keV gas temperatures and low (~ 0.2) best-fit abun- 
dances measured in NGC6861's tails. Long, hot tails can be 
formed by turbulent-viscous stripping of galaxy gas and its 
subsequent mixing and thermalization with the group IGM 
during supersonic motion of the galaxy through the ambient 
group gas, or by gravitational focusing of the IGM into an 
adiabatically heated Bondi-Hoyle wake (see, e.g. Machacek 
et al. 2005b). However, when such supersonic inf all o ccurs 
in the plane of the sky, as we expect here (see 33.31 ), we 
should see a rise in temperature and density along a nar- 
row shock front preceding the motion, followed by a cold 
front edge closer to the galaxy. We do see higher tempera- 
ture gas (kT ~ 1.3^02 keV), compared to the ambient IGM 
(~ 0.7-0.9keV), in the sheath surrounding NGC 6861 (re- 
gion FA in Fig. |T0b . If the factor ~ 1 .9 temperature rise in the 
sheath is due to a shock, we would expect to see a correspond- 
ing factor ~ 2 increase in density from the Rankine-Hugoniot 
conditions and, thus, a surface brightness increase by a factor 
4 in the same region (see, e.g. Landau & Lifshitz 1959). Al- 
though we find no evidence in these data for the sharp change 
in surface brightness (and thus density) that would be associ- 
ated with such a shock, nor do we see the associated merger 
cold front expected if NGC 6861 were moving supersonically 
in the plane of the sky through the ambient IGM, projection 
effects in these sparse data may make such features difficult 
to observe. 

A more promising scenario may be that surface brightness 
and temperature asymmetries are the result of a gravitational 
encounter between NGC 6868 and NGC 6861, each the dom- 
inant galaxy in their own subgroup. Then the dark matter dis- 
tribution is bimodal, with each galaxy residing near its own 
local minimum of the gravitational potential. This would ex- 
plain why double /3-models provide a good description of the 
mean surface brightness profile around each galaxy. However, 
since the merger is likely occuring in the plane of the sky, this 
subgroup structure may be difficult to discern from the opti- 
cal galaxy distributions. We may be viewing the encounter 
after the subgroups have passed through pericenter, but be- 



fore the gas halos have been completely stripped and merged. 
The strong temperature and surface brightness asymmetries 
around each of the merging partners may be due to non- 
hydostatic motions of the gas ('sloshing'), as gas is displaced 
from the center of each subcluster's dark matter potential, 
similar to those produced in simulations of premerger encoun- 
ters of more massive subcluster halos. When gas is present in 
both halos, cold fronts with cool spiral-shaped tails may form 
about one subcluster, while higher temperature shock-heated 
features may form around the other (see Figure 9 of Ascasi- 
bar & Markevitch 2006). These simulations suggest that the 
asymmetries persist for gigayears after the encounter. Thus 
the observation of sloshing signatures around each galaxy, 
even though they are separated by 207 kpc (sound crossing 
time of 500 Myr), is reasonable. 

4.3. The Case for AGN Feedback 

AGN activity and black hole growth may be influenced by 
either ram pressure, through backflow of partially stripped 
material onto the central regions of the galaxy, or during 
a close gravitational encounter between two galaxies, that 
would also produce sloshing. Since both NGC 6868 and 
NGC 6861 currently host weak radio sources at their nuclei 
consistent with AGN, we look for signatures in the X-ray gas 
that might signal past episodes of more energetic AGN activ- 
ity. 

4.3. 1 . 'Ghost' Cavities Near NGC 6868? 

Figures |2] and [TT] show a roughly circular deficit in X-ray 
surface brightness centered 4.'9 (in projection) to the north- 
west of NGC6868's nucleus. In the upper panel of Fig- 
ure [TT] we use 8 annular sectors, centered on the nucleus of 
NGC 6868 and with inner and outer radii of 168" (22.5 kpc) 
and 387"6 (5 1.9 kpc), respectively, to investigate more quan- 
titatively this northwest cavity. We denote the position of each 
sector by the mean angle 4> for the sector measured counter- 
clockwise from the +x axis (west). In the lower panel we 
plot the number of net source counts in each annular sector 
in the 0.5-2keV energy band as a function of sector posi- 
tion. Error bars represent la uncertainties in the net source 
counts. We find that the number of net counts (187 ± 19) in 
the sector coincident with the northwest cavity (<f> = 37.5°) 
is lower by > 3er than the counts in its neighboring sectors, 
(281.3 ±21.4 counts at (j) = 81.5° and 440.2 ±24.7 counts at 
= -8.5°[351.5 ]). 

The cavity candidate is a 'ghost' cavity, since no radio emis- 
sion has as yet been observed in this region, has low contrast 
and no well defined rim, similar to the radio-faint cavity in 
2A 0335+096 (Birzan et al. 2004). Features of this kind need 
to be confirmed in deeper X-ray observations and explored 
with lower frequency radio observations to detect the radio 
plasma expected to reside in their interiors. Furthermore, the 
projected spherical morphology is so mew hat puzzling given 
the evidence for gas motions in 34. II and 34. 21 that might dis- 
tort the cavity's shape. However, assuming the cavity is an 
evolved buoyant bubble filled with relativistic radio plasma 
from past AGN activity, we may use the geometry of the cav- 
ity and the mean pressure in the surrounding IGM to estimate 
the enthalpy (4pV) carried in the bubble. We calculate the 
mean pressure of gas that was evacuated to form the bubble 
from the density model in eq. [3] and the temperature of the 
IGM given in Table [4] Assuming a spherical bubble of radius 
11.7 kpc (87") oriented close to the plane of the sky, we find 
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FIG. 11. — (upper) 0.5 — 2keV blank-sky background subtracted, expo- 
sure corrected Chandra image of NGC 6868 with a panda region overlaid 
to identify X-ray surface brightness deficits (cavities) in the group IGM. 
1 pixel = 1"968 X 1 "968 and the image has been smoothed with a 6" 
Gaussian kernel. The annular sectors are centered on NGC6868's nucleus 
(20''09" , 54.1 J ,-48°22'46.4") with inner and outer radii of 168" (22.5 kpc) 
and 387" 6 (5 1 .9 kpc), respectively. North is up and east is to the left, (lower) 
Net source counts as a function of the mean angle cj> of the annular sector, 
measured counterclockwise from west. Error bars represent 1 a uncertainties 
in the source counts. Note the surface brightness deficit at <j> = 36.5°. 



the enthalpy carried by the bubble of ApV ~ 5.9 x 10 56 ergs, 
similar to ghost cavities in Abell 262 (Blanton et al. 2004) or 
HCG 62 (Birzan et al. 2004) and to younger cavities in galax- 
ies, such as M87 (Young et al. 2002; Forman et al. 2005) and 
NGC 4472 (Biller et al. 2004). 

The age of the bubble may be estimated by assuming that 
(1) the bubble rises with the speed of sound c s in the ambient 
medium t = d/c s , (2) the bubble rises buoyantly at the termi- 
nal velocity v, such that t =d/vt ~ dy/SC/2gV, where V is 
the volume of the bubble, S its cross-section, and C = 0.75 
the drag coefficient, or (3) the time t = 2^J rjg for gas to re- 
fill the displaced volume as the bubble moves upward, where 
r is the radius of the cavity and g the acceleration of grav- 
ity (Churazov et al. 2001; McNamara et al. 2000; Nulsen 
et al. 2002; Birzan et al. 2004). For the northwestern bub- 
ble these imply bubble ages of 88, 107, 119Myr, respec- 
tively, for the three methods, and a mechanical luminosity 
of Lmech ~ 1.8 x 10 41 ergs" 1 . Bubble ages of ~ 10 8 yr are 



typical for ghost bubbles in other systems (see, e.g. Birzan 
et al. 2004). 

A partially rimmed, spherical X-ray surface brightness fea- 
ture, with radius ~ 60" may also be seen in Figure [2] located 
to the southwest of NGC6868's nucleus. A filament of 
hot gas (region F in Fig. ITOb is seen to cross the candidate 
cavity. If this feature is also a buoyant ghost bubble, perhaps 
broken by the hydrodynamic motions in the gas caused by 
interactions between NGC 6868 and NGC 6861, we estimate 
its enthalpy, age and mechanical power to be 4 x 10 56 ergs, 
55 — 78Myr, and ~2x 10 41 ergs" 1 , respectively, comparable 
to that of the northwest cavity. 

4.3.2. AGN Heating in NGC 6861? 

The surface brightness in the plume rising to the northeast 
from NGC 6861 and in the western tail increase by factors ~ 2 
and ~ 3, compared to regions to the southeast at compara- 
ble distance from NGC 686 l's nucleus, and have a ~ 0.5 keV 
higher temperature than the central region of that galaxy (see 
Fig. [10] and regions FA and WT compared to G in Table [9]). 
One possible explanation is that gas in these regions may have 
been disturbed and heated by AGN activity. Assuming that 
the western tail is contained in a cylindrical volume oriented 
in the plane of the sky with radius 18 kpc and length 37.7 kpc 
(based on the WT spectral region), we find the tail contains 
~ 2.6 x 1O 9 ?7 1 ' 2 M0 of hot gas. Since the gas in the tail ap- 
pears clumpy, the filling factor -q is likely less than one, so 
that this is an upper bound on the actual mass of hot gas in 
the tail. The energy necessary to raise the temperature of the 
gas in the tail from ~ 0.7 keV found for gas in the bright cen- 
tral region of the galaxy (region G) or the surrounding IGM 
to 1.18kev measured in the tail is ~ 6 x 10 57 ergs. Outburst 
energies of up to 10 61 ergs have been measured from shocks 
observed in, for example, Hercules A (Nulsen et al. 2005a) 
and Hydra A (Nulsen et al. 2005b), while outburst energies 
~ 10 58 ergs have been measured in M87 (Forman et al. 2005), 
whose 3 x 10 9 M Q central black hole has mass comparable to 
that inferred for NGC 6861 from the Mbh - relation. Thus 
energetic outbursts from a 2.5 x 1O 9 M black hole at the cen- 
ter of NGC 6861 would have sufficient energy to heat and pro- 
duce turbulent motions in the gas, that might also contribute 
to the disturbed gas morphology observed near NGC 6861. 

5. ON NGC 6861'S BLACK HOLE MASS 

The discrepancy between the black hole mass prediction 
from the Magorrian relation (eq. O and the Mbh _ o* rela- 
tion (eq. [B is still difficult to understand. If NGC 686 l's 
black hole mass is large, then either an existing massive stel- 
lar bulge has been stripped from NGC 686 1 or the coevolution 
of NGC 686 1 's stellar bulge has been inhibited or delayed dur- 
ing a period of rapid black hole growth. In the first case, hy- 
drodynamic processes such as ram pressure and turbulent vis- 
cous stripping, plasma outflows from AGN, or gas sloshing 
may heat and redistribute the gas in and around NGC 6861, 
but would be ineffective at changing the galaxy's stellar dis- 
tribution. Stars in the outer halo could be stripped tidally 
if the impact parameter of the collision between NGC 6861 
with NGC 6868 is sufficiently small, as observed in dumbbell 
galaxy systems such as NGC 4782/NGC 4783. However, even 
in the deeply penetrating collisions of such dumbbell galax- 
ies, simulations show that only 10-20% of the stellar mass 
is expected to be stripped (Madejsky & Bien 1993). For a 
2.5 x 10 9 M Q black hole, equation [2] predicts a stellar bulge 
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mass for NGC 6861 of ~ 10 M©, such that previous interac- 
tions would have had to strip ~ 90% of the stellar mass to be 
in agreement with the K-band mass estimates given in TableQ] 
One would expect to see evidence for such cataclysmic tidal 
disruption in tidal bridges, tails or other stellar morphological 
distortions, that are not observed in the stellar distributions of 
either NGC 6861 or NGC 6868. Also this scenario would im- 
ply the direct passage of the subgroup gas cores through each 
other and much more disruption in the X-ray gas in NGC 6868 
than is observed. 

The second case, in which the growth of the stellar bulge 
is inhibited, may be more likely. NGC 6861 shows evidence 
for a stellar bar (Koprolin & Zeilinger 2000). Stellar bars are 
likely to be formed in mergers and allow angular momentum 
of infalling matter to be dissipated. Gas and dust may then 
be transported along the bar to the galaxy's nuclear region to 
fuel one or several episodes of rapid black hole growth. If 
this accretion event triggers an AGN outburst in NGC 6861, 
gas otherwise available for star formation could be heated 
and inflated, delaying the expected growth of the bulge (Silk 
& Rees 1998). The total energy deposited into the ambient 
gas from an episode of rapid black hole growth is feAMc 2 , 
where AMc 2 is the rest mass energy of the accreted matter, 
e is the fraction of the rest mass energy released in the ac- 
cretion event, and /, the feedback efficiency, measures how 
efficiently that energy is coupled to the surrounding galaxy 
and subgroup gas. For the standard assumption of radiative 
efficient accretion (e = 0.1) and efficient feedback / > 0.05 
(Di Matteo et al. 2008; Churazov et al. 2002), the total en- 
ergy deposited in the surrounding gas as the black hole grows 
from 2.1 x 1O 8 M predicted by the Margorrian relation to 
2.5 x 1O 9 M0, consistent with the Mbh-c* relation, would 
be large (>2x 10 61 ergs). This is comparable to the most en- 
ergetic AGN outbursts observed in massive clusters (see, e.g. 
Nulsen et al. 2005a for Hercules A; McNamara et al. 2002 
for MS 0735.6+7421; Nulsen et al. 2005b for Hydra A). If 
accretion proceeds at the Eddington rate, 

M Edd = (O.22/e)(M BH /lO 8 )M yr" 1 , (5) 

the accretion timescale would be ~ 1 10 Myr and the mean ra- 
diative power, ~ 10 47 erg s" 1 . The power that couples to the 
surrounding gas (~ 5 x 10 45 ergs _I ) is more than three or- 
ders of magnitude larger than that needed to balance X-ray 
cooling in the galaxy. Such an event would likely completely 
disrupt the gas in NGC 6861 and its subgroup. However the 
accretion and feedback efficiencies are highly uncertain. If the 
growth of NGC 686 l's black hole was triggered by the ongo- 
ing merger with the NGC 6868 subgroup in AS0851, either 
the accretion is radiatively inefficient (e << 0.1), the released 
rest-mass energy is only weakly coupled to the galaxy and 
subgroup gas (f << 0.05; see, e.g., Kurosawa et al. 2009), 
and/or the black hole grows through multiple, less energetic 
outburst events over a much longer timescale. For example, 
if ft ~ 10~ 6 , the energy from AGN activity deposited in the 
surrounding subgroup gas from this episode of rapid black 
hole growth would be ~ 10 58 ergs, sufficient to heat the west- 
ern tail. Alternatively, the period of rapid black hole growth 
may have occurred at very early times, before NGC 6861 was 
assembled into the current galaxy group/subgroup structures 
in AS0851. Then a very strong AGN outburst, as indicated 
above, may have swept gas from the galaxy, delaying bulge 
formation by quenching star formation and heating or blowing 



away any surrounding subgroup gas. Growth of the associated 
stellar bulge for NGC 6861 may only now be occurring, not 
through star formation, but by the cannibalizaton of the stel- 
lar distributions already in place as NGC 6861 merges with 
nearby galaxies in the group and ultimately NGC 6868 (Silk 
& Rees 1998). 

A third possibility is that the Magorrian relation holds and 
the central black hole in NGC 6861 is — 2 x 10 8 M Q , typi- 
cal for early type galaxies in groups. Then gravitational in- 
teractions between NGC 6861 and NGC 6868 or other group 
galaxies would also need to increase the central stellar ve- 
locity dispersion in NGC 6861 to anomalously high values. 
Strong tidal forces in a deeply penetrating encounter may dis- 
tort stellar velocity profiles, causing them to increase with 
distance from the galaxy center (see, e.g. Madejsky & Bien 
1993 for NGC 4782). For finite apertures, this might contam- 
inate the measurement of the central stellar velocity disper- 
sion and thus cause the black hole mass to be overestimated. 
For example, the measurement of 382 ±7kms _1 by Kopro- 
lin & Zeilinger (2000) is the average stellar velocity disper- 
sion within 5" (670pc) of the galactic center, well outside the 
sphere of influence of the supermassive black hole. However, 
if the velocity dispersion profile has been affected by tidal 
forces in a merger, as in NGC 4782, one would expect that the 
stellar velocity dispersion would increase when averaged over 
larger radii, which is not observed in NGC 6861 (Koprolin & 
Zeilinger 2000). It is also difficult to understand in this sce- 
nario why the stellar velocity dispersion of NGC 6861 would 
be so profoundly affected by the interaction, while that of its 
likely interacting partner NGC 6868 is not. If the Magorrian 
relation holds, the anomalously high stellar velocity disper- 
sion in NGC 6861 may require NGC 6861 to have experienced 
one or more recent minor mergers unrelated to the interac- 
tion with NGC 6868 that is responsible for the large scale gas 
sloshing observed in the group. 

6. CONCLUSIONS 

In this paper we used Chandra observations of NGC 6868 
and NGC 6861 in the AS0851 galaxy group to measure the 
density and temperature structure of gas in these two dom- 
inant group elliptical galaxies. This system is particularly 
important because, although the brightest group galaxy has 
a black hole mass M B h ~ 3.3 x 10 8 M Q predicted by the 
A/bh - cr* relation, that is both consistent with the Margor- 
rian MBH-A^buige relation and typical of black hole masses 
inferred for the dominant galaxies in other cool groups, the 
central black hole mass inferred from the Mbh - c* relation 
for the second brightest group galaxy NGC 6861 (Mbh ~ 
2.5 x 1O 9 M ) is more than an order of magnitude higher than 
that predicted by the Magorrian relation. The goal of this 
work was to use these X-ray data to explore possible inter- 
action scenarios (infall and gas stripping, gas sloshing, AGN 
outburst activity) that might affect galaxy evolution and black 
hole growth. 

We find: 

• The X-ray surface brightness distributions around 
NGC 6868 and NGC 6861 are each well described 
by double /3-models with inner (outer) beta and core 
radii of 0.65 and 0.49 kpc (0.52 and llkpc) for 
NGC6868 and 0.61 and 0.5 kpc (0.38 and 11. 9 kpc) 
for NGC 6861. This suggests that NGC 6868 and 
NGC 6861 may each be the dominant galaxy in a 
galaxy subgroup that is merging. Lack of evidence for 
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substructure in the galaxy reshift distribution suggests 
that we would be viewing such a subgroup merger close 
to the plane of the sky. 

X-ray surface brightness and temperature maps of 
NGC 6868 show a cold front 23 ± 3 kpc to the north 
of the galaxy nucleus and a spiral-shaped tail of cool 
0.62 keV gas to the south, suggesting relative motion 
between the galaxy and the surrounding IGM. Analy- 
sis of the density and temperature across the edge con- 
strains the relative motion between NGC 6868 and the 
IGM to be at most transonic (< 460 km s" 1 ). While 
we cannot rule out the possibility that the tail is com- 
posed of ram-pressure stripped galaxy gas, the tightly 
wound spiral-like morphology of the tail in the temper- 
ature map strongly suggests that that the tail is a result 
of the 'sloshing' of gas in the gravitational potential, 
set in motion by the previous gravitational encounter 
between the NGC 6868 and NGC 6861 subgroups. 

X-ray surface brightness and temperature maps for 
NGC 6861, 206 kpc to the west of NGC 6868, show 
a complex morphology also suggestive of a subgroup 
merger. While spectra of the central ~ 5 kpc of 
NGC 6861 containing the bright stellar disk is well de- 
scribed by a cool 0.66^'^ keV thermal plasma model 
plus a power law with photon index -1 .87 for the cen- 
tral active galactic nucleus (AGN), the temperature of 
gas in the 5- 14 kpc 'sheath' outside NGC6861's cen- 
tral region is a factor ~ 2 higher. A hot (0.9- 1.2keV), 
bifurcated tail is seen trailing NGC 6861 38 kpc to the 
west and northwest. These features may suggest that 
the NGC 6861 subgroup is undergoing ram-pressure 
or turbulent-viscous stripping during supersonic infall 
towards NGC 6868. However, comparing the den- 
sity and temperature morphologies of NGC 6868 and 
NGC 6861 taken together with simulations from the lit- 
erature, the more likely explanation is that the hot fea- 
tures in NGC 6861 are the consequence of gas sloshing 
caused by the ongoing merger of two subgroups, whose 
respective hot gas halos have not yet been stripped. 



• NGC 6868 may have undergone a period of merger- 
induced AGN activity. X-ray surface brightness im- 
ages show weak evidence for roughly circular surface 
brightness deficits to the northwest and southwest of 
NGC 6868. If confirmed in deeper X-ray images or 
low frequency radio observations to be buoyant bubbles 
filled with relativistic plasma from previous AGN activ- 
ity, the bubbles' enthalpies (5 x 10 56 ergs), mechanical 
power (2 x 10 41 erg s" 1 ) and ages (~ 60- lOOMyr) are 
typical for 'ghost' cavities observed in other systems. 

While the X-ray data favor a sloshing interaction sce- 
nario aided by AGN activity to explain possible cavities near 
NGC 6868 and gas heating in/near NGC 6861, an explana- 
tion for the discrepancy between the Magorrian and Mbh - 
relation predictions for the black hole mass at the center of 
NGC 6861 remains challenging. Deeper X-ray observations 
that probe closer to the central engine and map with greater 
precision the complex density and temperature structure of 
the gas are needed to constrain the gasdynamics and orbital 
parameters of the merger and determine the energies and 
timescales of associated AGN activity. More detailed theo- 
retical modeling in gravity + hydrodynamic numerical simu- 
lations that are guided by those data may then be able to deter- 
mine how mergers affect the high mass end of the black hole 
mass relations and solve the mystery of NGC 6861 's black 
hole mass. 
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TABLE 1 

Galaxy Properties and Black Hole Masses 



Name 


^ (total) 








M BH 






(10"M Q ) 


(kms" 1 ) 


(1O 8 M ) 


(1O 8 M ) 


NGC 6868 


7.317 


1.83 


250 ±10 


3.3 


3.1 


NGC 6861 


7.708 


1.28 


414 ±17 


25 


2.1 



NOTE. — Columns (1) galaxy name; (2) total 2MASS K-band luminosity (NED); (3) stellar mass (bulge mass) estimated using K-band luminosities 
(Gilfanov 2004); (4) stellar velocity dispersion (Wegner et al. 2003); (5) and (6) central black hole mass from the Mbh — (eq. [B Tremaine et al. 2002) and 
the MBH - Mb u i e e (eq. [2] Haring & Rix 2004) relations, respectively. 



TABLE 2 

Density Model Fits to NGC 6868's Edge 



Fit range 
(kpc) 


"1 


' edse 

(kpc) 


J 


"I 

(10~ 3 cnT 3 ) 


n 2 

(10~ 3 cnT 3 ) 


[7.9.34.0] 


-0.40 


22.7 


2.2 


1.52 


0.84 


[7.9,68.0] 


-0.15 


20.3 


1.8 


1.47 


1.00 


[5.8,34.0] 


-0.82 


26.5 


2.6 


1.38 


0.66 


[5.8,68.0] 


-0.64 


22.5 


1.7 


1.16 


0.85 


NOTE. — Density models fits across the 0.5 - 
«2 are electron densities inside and outside the ( 
the edge, respectively. See Table|4] 


2 keV surface brij 
;dge, respectively, 


;htness edge in NGC 6868 shown in Figs. [2]and[8]for the density model given in eq. f5] n\ and 
assuming (kT,A) of (0.64 keV, 0.5Zq) and (0.72 keV, 0.3Zq) inside and to the north outside 
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TABLE 3 

NGC 6868 Spectral Regions Across the Edge 



label 


shape 


center 


radius" 


angle 






(RA, Dec) 


(arcsec) 


(deg) 


NEdge_! 


elliptical sector 


20:09:54.1,-48 : 22 : 46.4 


82.4,131.9 


6 


NEdge+o 


elliptical sector 


20:09:54.1,-48 : 22 : 46.4 


131.9,211. 


6 


NEdge +1 


elliptical sector 


20:09:54.1,-48 : 22 : 46.4 


211. ,337.6 


6 


S+i 


elliptical sector 


20:09:54.1,-48 : 22 : 46.4 


211. ,337.6 


6 



NOTE. — NEdge regions are concentric to the bounding ellipse and lie within the sector shown in Fig. |7]and described in 34.11 S + i uses the same elliptical 
annulus as NEdge+i but is constrained to lie between 213° and 327.5° (to the south). All angles are measured counter-clockwise from west. Coordinates are 
J2000. " Radii are semi-minor inner (outer) radii for elliptical annular sectors. 



TABLE 4 

NGC 6868 Spectral Models Across the Edge 



region 


counts 
(keV) 


kT 

(Z©) 


A 


r 


X 2 /doi 


( A") 


NEdge_! 
NEdge_! 
NEdge_[ 
NEdge_i 
NEdge-n, 
NEd ge+<) 
NEdge +1 
NEdge +1 
S + i 


313 
313 
313 
313 
286 
286 
464 
464 
691 


0.65+os 
0.65 ± 0.07 
0.64 ±0.08 
0.66 ±0.06 

072+o.os 

73+$ ;| W 

-Q-QS 



u - oy -o.io 

63+ 05 
U - OJ -0.06 


0.3 

0.5 

0.5 

1.0 

0.2 

0.3 

0.3 
11+013 
29+ 24 


1.6 


11.6/19 
11.2/19 
10.5/18 
11.1/19 
30.2/31 
30.6/31 
65.7/55 
61.7/54 
72.1/68 


6.2 

9.2 

9.2 

16.8 

4.7 

6.1 

6.1 

3.3 

6.1 



NOTE. — Spectral models are an absorbed APEC model or absorbed APEC + power law model. Hydrogen absorption is fixed at the Galactic value 
(Nh = 3.9 X 10 20 cm~ 2 , Kalberla et al. 2005). Parameters without 90% CL uncertainties are fixed. "A is the X-ray emissivity in the 0.5 — 2 keV energy band in 
units of 10 -15 photon cm 3 s -1 . 



TABLE 5 
NGC 6868 Edge Analysis 



case" 


Ti/T 2 


(Ai/A 2 )'/ 2 


«l/«2 


Pl/P2 


Mach 


velocity 
(kms- 1 ) 


1 


0.89 


0.99 


1.67-2.56 


1.49-2.28 


0.72-1.08 


318-475 


2 


0.88 


0.81 


1.37-2.10 


1.21-1.85 


0.47-0.9 


207-397 


3 


0.90 


0.60 


1.01-1.55 


0.90-1.39 


0-0.65 


0-287 



NOTE. — Case (1,2,3) correspond to spectral models given in Table[4]with abundances inside the edge of (0.3,0.5,1.0), respectively, and temperature, 
abundance and sound speed outside the edge of 0.72 keV, 0.3 Zq and 440 km s . 



TABLE 6 

NGC 6868 Spectral Regions for Temperature Map Features 



label 


shape 






center 




radius" 


angle 










(RA, Dec) 




(arcsec) 


(deg) 


jb 


ellipse 


20 


10 


13.3,-48 : 27 


16.1 


196, 105 


333 




ellipse 


20 


09 


40.0,-48 :27 


02.4 


149.6,92.8 


21.5 


F 


ellipse 


20 


09 


45.4,-48 : 24 


27.0 


150., 62.7 


29.2 


C 


circle 


20 


10 


04.9,-48 : 22 


54.4 


72.5 




N 


circle 


20 


09 


53.7,-48 : 22 


46.6 


6.2 





NOTE. — All angles are measured counter-clockwise from west. Coordinates are J2000. " radii are semi-(major,minor) axes for ellipses and radius for circles 
* T is the sum of the two elliptical regions given. 



Merger of NGC 6868 and NGC 6861 
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TABLE 7 



NGC 6868 Spectral Models for Temperature Map Features 


region 


counts 


kT 


A 


r 


X 2 /dof 




(keV) 


(Z©) 








T 


785 


0.62 ±0.04 


4+1-07 
u - M).20 




79.0/79 


F 


607 


0.71 ±0.05 


0.3 




37.3/37 


F 


607 


0.70 ±0.05 


0.2 




35.2/37 


C 


504 




0.5 




38.2/28 


N 


790 


63+ - 07 


1.0 


1.6 


37.7/31 


N 


790 


o 63+ - 07 

-Q-Bfi 


1.0 


1.4 


32.1/31 


N 


790 




1.0 


2.0 


51.9/31 



NOTE. — Spectral models are an absorbed APEC model or absorbed APEC + power law. Hydrogen absorption is fixed at the Galactic value Nh = 3.9 X 
10 20 cnT 2 , (Kalberla et al. 2005). Parameters without 90% CL uncertainties are fixed. 



TABLE 8 
NGC 6861 Spectral Regions 



label 


shape 


center 


radius" 


angle 




(RA, Dec) 


(arcsec) 


(deg) 


E 


rectangle 


20:07:51.8.-48 : 22: 01.0 


284.7,392.3 





G 


ellipse 


20:07: 19.5.-48 : 22 : 12.9 


31.6,19.4 


33.8 


FA* 


annular sector 


20 :07 : 19.5.-48 : 22: 12.9 


39.9,106.8 




WT 


rectangle 


20:07:01.6,-48:22: 13.7 


281.4,269.6 





NW 


rectangle 


20:07:14.3,^18:17:46.5 


372.8, 147.2 


47.3 



NOTE. — Angles are measured counter-clockwise from west. Coordinates are J2000. " radii are (length, width) for rectangles, semi-(major,minor) axes for 
the ellipse, and (inner.outer) radius for the annular sector. * annular region is constrained to lie between 80° and 282° measured counterclockwise from west. 



TABLE 9 
NGC 6861 Spectral Models 



region 


source 


kT 


A 


r 


x 2 /dof 




counts 


(keV) 


(Z©) 







El 

E2 

G 
FA 
WT 
NW 



404 
404 

436 
214 
621 

385 



61+0- 24 

-0.07 

1 26+0- 44 

J - 4 °-Q.26 
52+o.n 

066+ffl 

1.1 8 |]« 
0.90 ±0.09 



027 +<) - 04 

0.2 

0.2 

0.5 

5 +L8 

-0.28, 
25+ 421 
-013 
n 72+0-26 



28.5/22 

26.3/21 
8.6/16 
10.6/9 

35.7/42 
39/28 



NOTE. — Spectral models are an absorbed APEC thermal plasma model, two component absorbed APEC or absorbed APEC + power law. Hydrogen 
absorption is fixed at the Galactic value Nh = 3.9 X 10 20 cm" 2 , (Kalberla et al. 2005). Parameters without 90% CL uncertainties are fixed. 



